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Abstract 

In this paper we explore the possibihty of a steady-state entanglement of two two-level atoms 
inside a pumped cavity by taking into account cavity leakage and the spontaneous emission of 
photons by the atoms. We describe the system in the dressed state picture in which the coherence 
is built into the dressed states while transitions between the dressed states are incoherent. Our 
model assumes the vacuum Rabi splitting of the dressed states to be much larger than any of 
the decay parameters of the system which allows atom-field coherence to build up before any 
decay process takes over. We show that, under our model, a pumping field cannot entangle two 
closed two-level atoms inside the cavity in the steady-state, but a steady-state entanglement can 
be achieved with two open two-level atoms. 
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I. INTRODUCTION 



In this paper we investigate the steady-state entanglement of two two-level atoms inside 
a high-Q cavity which is pumped by an external field. We assume a large vacuum Rabi 
splitting, and take into account the presence of cavity leakage and spontaneous emission by 
the atoms. We demonstrate that a pumping field cannot entangle two two-level atoms inside 
a cavity in the steady-state, but it is possible to have steady state entanglement between 
two open two-level atoms in a cavity in the presence of a pumping field. This is achieved by 
optically pumping half of the atomic population outside of the decoherence-free subspace of 
the system, and the process yields an entanglement equivalent to a Bell state content of |. 

Entanglement between two partitions of a system is often due to some constraint placed 
on the dynamics of the system. This constraint can take the form of energy conservation, 
momentum conservation, or structural constraint to name a few. In this paper we direct 
our attention to the entanglement that arises from the constraints on the energy, or more 
specifically, the number of excitations in our system. 

In a simple system of two two-level atoms inside a cavity, initially in their ground states, 
and with one photon in the cavity mode, we know that the two atoms get entangled and 
disentangled periodically in time. This could lead one to think that if the two atoms were 
placed in a sufficiently high-Q cavity with small losses introduced and externally pumped 
such that the average photon number in the cavity is less than or equal to one, then the two 
atoms would still get entangled in some fashion. In the simplest case in which one pumps 
the cavity with a single photon source, and constrains the pumping rate to be less than or 
equal to the decay rate of the system, it would seem plausible that the two atoms inside 
the cavity can be entangled. In fact, making a weak field assumption and only taking into 
account the first order correction to the reduced density matrix of the atoms will yield a 
non-zero entanglement between the atoms. However, we show that if one were to take into 
account the higher order corrections in the density matrix, then the entanglement between 
the atoms vanish. The result is due to the fact that the concurrence measure is a nonlinear 
function of the parameters of the atomic density matrix. The weak field assumption is only 
valid if one is calculating some property of the atoms which is linearly dependent on the 
parameters. 

Recent experimental advances in atomic traps and cavity QED 
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offer exciting 



2 



new possibilities in quantum computation and quantum networking using the entanglement 
between atoms in the cavity . It is within our technological limits to trap and manipulate 

n 

atoms inside a microcavity to study their entanglement behavior [3]. Therefore, it is impor- 
tant to characterize the entanglement of the atoms in the cavity, and how the entanglement 
is changed by the dynamics of the system. 

Theoretical work on atoms inside a cavity , [lO| suggests not only entanglement 



between atoms in a cavity, but the possibility of a steady-state atomic entanglement in 
presence of a weak pump. Here we will investigate under what circumstances the atoms 
inside the cavity will get entangled with each other, and whether a steady-state entanglement 
can be attained. 

Our approach to characterizing the entanglement between atoms inside a cavity is through 
the dressed state picture of the atom-field states in the limit of large vacuum Rabi splitting. 
This approach enables us to use rate equations to characterize the transitions between the 



dressed states 
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121]. All the coherent effects are contained within the dressed states. 



and the transitions between the dressed levels are incoherent effects which are attributed 
to spontaneous decay of the atom, cavity leakage, or pumping. Also, for simplicity, we 
ignore longitudinal dipole-dipole interaction between the atoms and motional effects in our 
calculation. The advantage to our approach, which expresses the density matrix of the 
system as a mixture of the dressed states, is that we can identify which dressed states 
contribute to the entanglement between the atoms. This will more directly suggest how one 
may manipulate the system in order to maximize entanglement between the atoms. We will 
show that in the steady-state regime a pumping field cannot entangle two closed two-level 
atoms in the steady-state, but it is possible to entangle two open two-level atoms in the 
steady-state with an entanglement content worth one half a Bell state through optically 
pumping atomic population out of the decoherence-free subspace of the system using a 
pumping field. 

First, we will investigate the case in which two closed two- level atoms are inside a cavity 
by taking into account the n = 0, n = 1, and n = 2 dressed states, and we show there is no 
entanglement between the atoms. Then we will show that a calculation involving only the 
n = and n = 1 dressed states yield a non-zero entanglement between the atoms. Although 
this approach is reasonable for weak pumping fields when one is seeking the density matrix for 
the purpose of determining the population distribution, or other characteristics of the atoms 
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which depend hnearly on the density matrix, we will show that taking into account only the 
n = 1 dressed states is insufficient for the purposes of determining the entanglement between 
the atoms. We then proceed to show that a steady-state entanglement can be obtained in 
the presence of a pumping field for the case of two open two-level atoms even when we 
consider all the dressed states of the system. 

II. MODEL SYSTEM: TWO CLOSED TWO-LEVEL ATOMS INSIDE A CAVITY 
A. Multiple Excitation in the Cavity 

In this section we will show that a pumping field cannot generate entanglement between 
two two-level atoms inside a high-Q cavity in the presence of cavity leakage and spontaneous 
emission by the atoms. The atoms are both initially in the ground state, and they are 
placed in the cavity with the cavity mode being on resonance with the atomic transition. 
The external pumping field is also on resonance with the atomic transition in our model. 
Furthermore, to be consistent with the high-Q cavity assumption, we will stipulate that the 
Rabi splitting is much larger than the spontaneous decay rate and the cavity leakage rate. 
The large Rabi splitting justifies the rate equation approach we will take to analyze our 
system. 

The system we are considering consists of two atoms, each with energy structure shown 
in Fig.([T]), in a cavity which is externally pumped on resonance with the atomic transition. 
The system can lose energy through cavity leakage or through spontaneous emission of the 
atoms. 




FIG. 1: Two closed two- level atoms inside a cavity 
Since we wish to characterize this system using rate equations, we will ffist obtain the 
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dressed states of the closed system, and later calculate the transition rates between these 
dressed states. We assume that the Rabi splitting in this system is large such that the 
coherence of the atoms and the field is contained in the dressed states. Therefore, pumping 
population from the ground state up to these dressed states results in populating an inco- 
herent mixture of these dressed states. We will study the entanglement between the atoms 
in the steady-state regime after we obtain the rate equations that govern the population 
distribution of the system. 

The Hamiltonian of the closed system in the interaction picture is 

Hi = g^ai^a^ + ^s^g^a^ + H.C. (1) 

where a^^^ is the atomic transition operator for the n*'* atom and is the field creation 
operator. For the sake of simplicity we will assume the coupling constants, gi, to be the 
same for each atom. In the case in which there is only one excitation in the system there 
are three essential states, 

|11;1) 

|12;0) (2) 
|21;0), 

and three dressed states, 

IXo) = ;^(|12;0)-|21;0)), (3) 
|X+) = ^|11;1) + ^(|12;0) + |21;0)), (4) 

|X-) = ^|11;1)-^(|12;0) + |21;0)). (5) 

Here \ab;c) — |a)i ig) |6)2 <E) |c)/ indicates the first atom is in state |a), the second atom in 
state 1 6), and the field in state |c). 

In order to construct the rate equations governing the distribution of the populations 
we will now need to calculate the spontaneous emission rate, pumping rate, and the cavity 
leakage rate of the system. The spontaneous emission rate between states \a) and \b) is 
proportional to |(a|(i(j6|&)P where dab is the dipole matrix element between states \a) and 
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\b). Therefore, the rate of decay of the dressed states into the ground state, \g) — |11; 0), is 
given by. 



7+^, = r|(x+|rf« + d(^)|^?)p = r, (6) 
7-., = r|(x-|rfW + rf(2)|^?)p = r, 

where F is the Einstein A coefficient of the 2 — >■ 1 transition of the single atom in free space. 

Likewise, the pumping rate of the cavity is proportional to |(c|d^|c— where |c) is 
the field state in the Fock basis. The pumping rates for the dressed states from the ground 
state is given by, 

7r,^+ = n|(x+|at|^)|2 = n (7) 

t|n\|2 = n 



7r,_ = n|(x-|at|(?)| 



2 ' 



where FI is the single photon pumping rate inside the cavity (i.e. FI = HsourceT'^ where T 
is the amplitude transmission coefficient of the input cavity mirror and Hgource is the single 
photon emission rate of the pumping source) . 

The cavity leakage rate is obtained in a similar fashion, 

Ko^g^ K\{xo\a''\g)\'^ ^0, 

K+^, = K|(x+|at|^)|2 = f, (8) 

K.^^ = K\{x-\a^\g)\' = §. 
Here K is the power transmission coefficient of the cavity output mirror. 

The n >2 excitation of the system will have a different set of dressed states since there 
are four essential states in this case. The four essential states for n > 2 are. 



11; n) 




12; n- 


1) 


21;n- 


1) 


22; n- 


2) 



(9) 



The dressed states are then given by. 



|C) = ^(|12;n-l)-|21;n-l)), (10) 
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|C,) = -^(|ll;n)-|22;n-2)), (11) 
\4>l) = ^(|ll;n) + |12;n - 1) + |21;n- 1) + |22;n- 2)), (12) 

\(tr) = ^(|ll;n) - 1 12; n- 1) - |21;n- 1) + |22;n- 2)). (13) 

After we calculate the pumping rate, cavity leakage rate, and spontaneous emission rate 
by the atoms between the dressed states, as we have done before for the n — 1 dressed 
states and the ground state, we can construct the rate equation governing the population 
distribution amongst the n = 0, n = 1, and n = 2 dressed states. The rate equation is given 

by, 

Psi = -rp,i - f p,i + uPg + |rp,2 + rp<,.,2 + k{p,2 + Po',2) + fnPsi, 
p,2 = -irp,2-i^p.2 + np,i, 
p„,,2 = -rPo^2 - KPo' + |np,i, 

^ (14) 

Po = 0, 

Po,2 = 0, 

P,1 = P+ + P_, 

Ps2 — P+,2 + P-,2, 

where Pg is the population of the ground state of the system, P± is the population in the 
\x±) dressed states, Pq is the population in the \xo) dressed state, P±^„ is the population in 
\(f>±), Po,n is the population in the and Po',n is the population in Icf)^/). Here we assume 
that there is no population initially in the \xo) and I^JJ) dressed states, and because these 
states don't couple to any other states, they will not accumulate any population at later 
times. 

The steady-state solution to these rate equations is. 
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Pg = Ar(3r + 2K){2V + K){V + K), 
Psi=M2Yi{T + K){'iT + 2K), 

P,2=MAIi\T + K), (15) 
=Arn2(3r + 2K), 

M = 7TIP + Qii^K + enr^ + wukt + auk^ + er^ + ukt^ + 9k^t + 2K^. 

Are the two atoms in the cavity entangled? In order to answer this question we will 
have to calculate the entanErlement content of this bipartite two-level mixed state of the 

n 

system. Here we will employ the concurrence measure put forward by Wootters pj[. The 
concurrence measure of entanglement is defined as, 



C = max(^/\^ — \/A2 — \/A3 — a/A4, 0) (16) 

where Aj are the eigenvalues, in descending order of value, of the matrix pp (p = {ay ® 
ay) p* {ay® ay)). 

To determine the entanglement between the atoms in the cavity we will trace out the 
field component of the density matrix to obtain the reduced density matrix of the system. 
The reduced density matrix of the two atoms is given by. 



Patoms 

= PgPg + Pslpsl + Ps2Ps2 + Po',2Po',2 (17) 

where 
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(18) 



The important point to note in the dressed states for n > 2 is that, aside from the 10") 
which does not couple to any of the other dressed states through spontaneous emission, 
pumping, or leakage, the other three dressed 1^+), and, have no entanglement (i.e 
C = 0) between the atoms. This is in contrast with the n = 1 dressed states, |x±), which 
does contain entanglement between the two atoms {C = \). This suggests that any strong 
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pumping which puts appreciable population in the n > 2 will disentangle the atoms. On the 
other hand, it may seem to suggest that in the limit of weak pumping, in which the n > 2 
dressed states are not appreciably populated, the atoms may turn out entangled. We will 
show that the atoms do not get entangled no matter how weak the pumping field. 
The concurrence of Patoms is, 

C atoms = max{^iPsl + Ps2) - l\/iPs2 + 2P„, 2) (2Po',2 + Ps2 + ^Pg + 2P.l), O). (19) 

It is not easy to determine by observation if there exists some combination of parameters 
which would provide a non-zero concurrence. The plot of the expressions inside the max 
function of the concurrence (which we will call C) against 11 and K (in units of F) Fig. ([2]) 
suggests that there is indeed no combination of parameters which will provide a non-zero 
entanglement content. 




FIG. 2: Plot of C against 11 and K (in units of F) 

Can there be a miniscule amount of entanglement that we don't see in the plot? In order 
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to answer this let us take a closer look at the expression for concurrence. First, we will 
determine if there are any extremums in the expression for concurrence by looking at the 
partial derivatives of C with respect to the Pi parameters and setting the expression to zero. 
After differentiating C we obtain, 

|| = -4(8Po,,2+4P,2)=0, 
fe = |-4(4i'o',2 + 2P,2)=0, 

^ = i - f (4P, + 2P,i + 2Ps2 + 4P„, 2) = 0, (20) 
^ = -4(4P,2 + 8Pg + 4P,i + 8P„,,2) = 0, 
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It is easy to see that ^ = and = cannot be satisfied by positive real values of the 
Pi. The Pi have to be positive real values since they are the probability weighting of the 
density matrix, and the steady-state solution suggests a non-zero value for all the Pj. The 
= ^iid = takes some algebra, but one can also show that the equations cannot 
be satisfied within the restrictions of the parameters. Therefore, there are no extrema 
or saddle points in the concurrence, which is consistent with Fig. ([2]). This means that 
the extremum value lies on the boundary of the domain of the Pj. The four boundary 
points are {Pg = l,P,i = 0, P,2 = 0, Po,,2 = 0), {Pg = 0,P,i = 1,P,2 = 0, Po,,2 = 0), 
(Pg = 0, Psi = 0, P,2 = 1, Po',2 = 0), and (Pg = 0, P,i = 0, P,2 = 0, Po',2 = 1), of which only 
the {Pg = 0, Psi = 1, Ps2 = 0, Po',2 = 0) point offers a nonzero concurrence of C = ^. Hence, 
the only possible non-zero value of the concurrence is along the P^i axis from {Pg = 0, P^i = 
l,Ps2 = 0, Po',2 = 0). Our task is to determine the maximum value of Pgi constrained by 
the steady-state solutions to the rate equations. 

Using numerical techniques one can determine the maximum value of Psi to be Pgi ~ 
0.366 which implies Pg ^ 0.317, P,2 ~ 0.211, and Po',2 ~ 0.106 (see Fig. ([3])). These 
population yields a concurrence of C = max (v/(0.083) - 2v/(0.0640),0) = 0. This means 
that a pumping field that takes any population beyond the n = 1 dressed state leaves the 
atoms completely unentangled. 

There are a couple points here worth noting. First, we only extended the rate equations 
to include the n = 2 dressed states, and yet the concurrence suddenly dropped to zero. In 
reahty, the n > 2 dressed states do get populated as well, so as far as the entanglement of 
the atoms are concerned we are actually overestimating the concurrence. This is because 
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FIG. 3: Plot of Psi against 11 and K (in units of F) 



in reality only some fraction of that population resides in the n = 2 states and the rest is 
somewhere in the n > 2 dressed states. Therefore, this implies even less population in the 
\x±) dressed states, and, as a consequence, less entanglement between the atoms. Second, 
we would expect that in the limit of a weak pump the n = 2 dressed states would not 
get appreciably populated, so the entanglement between the atoms ought to be non-zero. 
However, our analysis indicates that any population at all in the n = 2 dressed states 
disentangles the atoms. Lets take a numerical example to illustrate this peculiar point. For 
n = .447 and K = 1 (in units of F) we get Pg = 0.97337, P,i = 0.02595, P,2 = 0.00042, 
and Po',2 = 0.00026. One would be tempted to neglect the Ps2 and Po'2 terms in the density 
matrix since they are two orders of magnitude smaller than P^i. If we only keep the leading 
two terms then we get a concurrence of C = .01. However, if we were to keep all the 
terms in the density matrix we would get a concurrence of C = 0. Although Ps2 and Po'2 
were two orders of magnitude smaller, by keeping these terms in the density matrix we 
see there is no entanglement between the atoms. In principle, we can always extend the 
accuracy of our model by increasing the Rabi splitting of the dressed states. Therefore, the 
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miniscule population in n = 2 disentangling the atoms is indeed the case for the model we are 
considering. Because concurrence is not necessarily linear in the parameters of the density 
matrix, one cannot ignore terms even when they are several orders of magnitude smaller. In 
order to get a non-zero entanglement between the atoms we have to put more population in 
the \x±) dressed states, and minimize the population transfer to the n > 2 dressed states. 
The only realistic option we have is to have cavity leakage rate that is nonlinear in the photon 
number since F is a fixed parameter of the atoms and 11 is a pumping parameter that 

does not differentiate between the cavity field excitation number. If the mirror of the cavity 
demonstrates a strong enough nonlinearity between the one excitation and two excitation 
of the cavity field then one can obtain a non-zero concurrence between the two atoms. 

Finally, we would like to offer a more intuitive justification for our use of the rate equa- 
tions. Clearly, the way we construct the rate equations suggests that the populations of 
the dressed states depends only on the populations of the other dressed states and not the 
coherence between the dressed states. The downward transitions, spontaneous emission of 
the atoms and cavity decay, are incoherent processes, so they do not impart any coherence 
onto the system. The upward transition, the external pumping, also does not impart any 
coherence in the steady-state limit. When a photon impinges on an atom initially in the 
ground state inside a cavity, the system is put in a coherent superposition of the dressed 
states whose amphtudes oscillate with the associated Rabi frequency. However, when a 
decay process is introduced, the atom-cavity system can lose the photon and go back into 
the ground state. At this point they wait for the next pump photon to come and restart 
the Rabi oscillation. In effect, the quantum jump that is introduced randomizes the initial 
phase of the dressed states Rabi oscillation. We do not know when the photon was lost 
nor do we know when the new pump photon has arrived. Therefore, in the longtime limit, 
after many decay and pumping processes have taken place, we do not expect to have any 
coherence between the dressed states since we will have no information concerning the phase 
of the Rabi oscillation. For this reason we can use rate equations to describe our system. 

B. Single Excitation in the Cavity 

In the previous section we demonstrated that a pumping field cannot entangle the atoms 
in the steady-state inside the cavity under our model. In this section we will briefly exam 
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what one might obtain for the entanglement between the atoms if one employs the weak 
field assumption and decides to only consider the n = 1 dressed states. The model seems 
plausible in the limit of a weak pumping field, but we will show how the weak pumping 
assumption does not translate to only considering the n = 1 dressed states. 

If one were to consider only the n = 1 dressed states then the rate equation governing 
the population distribution given by. 



Po 



0, 



P-) + f{P 



P- = -rP- - fP- + ^Pg- 



Up 



(21) 



CO 



n/2 



K/2,r 



n/2 




K/2,r 



4- |11;0> 



FIG. 4: Schematic of the dressed state 
Provided there is no population in Po sX t = 0, the steady-state solution of these equations 



IS, 



P4 



p~ = 



F, = iV(f + r), 



(22) 



r+f+n 

The total density matrix (including the two atoms and the field) is given by. 



p = Pg\g){g\ + P+\x+){x+\ + P-\x-){X-\ + Po\Xo){Xo 



(23) 
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Since we are only interested in the entanglement between the two atoms, we will trace out 
the field part of the density matrix to obtain the reduced density matrix for the two atoms. 
This is given by, 

TTf,M = + \ip^ + + \{p^ + p-)\^^){i^^\ + Po\r){r\ (24) 

where IV'^) = ^(|12) + |21)) and |^-) = ^(|12)-|21)) are the two Bell states. Substituting 
our previous steady-state solution to the above equation yields, 

Patoms = + r)|ll)(ll| + ||^+)(^+|] . (25) 

The concurrence for this reduced density matrix is simply C — N^. This model suggests 
that whenever the pump is on there is always some entanglement between the atoms inside 
the cavity. This result is contrary to the one we have obtained earlier. If we only consider 
the n — 1 dressed states, then there is always some entanglement between the two atoms 
in the steady-state when the pump is on. However, we have shown earher that, no matter 
how weak the pumping field, the two atoms do not get entangled in the steady-state. The 
discrepancy between the two results is due to the fact that concurrence is not a linear function 
of the density matrix. If concurrence were linearly related to the density matrix, then the 
weak field assumption would be valid, and the two atoms would indeed be entangled. But 
because concurrence is not a linear function of the density matrix (or, the wave function 
of the system), one cannot dismiss the n > 2 terms in the density matrix even when the 
population in those states are significantly smaller than the n = 1 and n = dressed states. 
Therefore, the weak field assumption does not translate to just keeping the lowest order 
terms in the density matrix if one is interested in calculating the entanglement content of a 
system. 

In our analysis we have assumed that the atoms were on resonance with the cavity field 
mode and the coupling constants of the atoms were the same. If the atoms had a different 
coupling constant gi, or detuned from the cavity field by Aj, this would result in the reduction 
of the entanglement between the atoms. The entanglement between the atoms here is due 
to the l'^"'") in the |x±) dressed states. Different coupling constants or detuning would turn 
the |^+) into --^^====(q;|12) -|-/3|21)) in the dressed states, and in general, (3. Because 
maximum entanglement is attained when — P — the effect of differing or Aj's 
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would reduce the entanglement between the atoms. Therefore, the case we are considering 
is the case for maximum possible entanglement between the atoms. 



III. MODEL SYSTEM: TWO OPEN TWO-LEVEL ATOMS INSIDE A CAVITY 

So, is there a way to generate steady-state entanglement with pumping/leakage through 
the cavity and spontaneous emission? Indeed there is a way to generate entanglement under 
these conditions using two open two- level atoms. The system we consider. Fig. ([5]), is two 
open two-level atoms inside a cavity with spontaneous emission down to both ground states 
and the cavity being in resonance with only one transition of the open two-level atom. 




FIG. 5: Two open two- level atoms inside a cavity 

Since the pumping field does not interact with |3), the Hamiltonian of the closed system 
is the same as the two closed two-level system model, and therefore the dressed states are 
also the same. The only difference is that when we construct our rate equations we have the 
possibility of the atom to fall into the |3) state through spontaneous emission. 

Now, if we were to start out in the |11;0) state as we did before then we will get no 
entanglement between the atoms since the |33; n) state is clearly separable, and we saw in the 
previous section that the population distribution in the dressed states yields a concurrence 
of zero. Here we will take advantage of the maximally entangled \xo) and states by 
setting the initial state of the two atoms to be a linear combination of the symmetric and 
antisymmetric states. Assuming we start out in the |11; 0) state, we will apply a vr-pulse to 
one of the atoms before we turn the pumping field on. This will take the |11; 0) state to the 
1 12; 0) state. In terms of the dressed states this is expressed as, 

|12;0) = ^(|x+)-|X-)) + ^IXo). (26) 
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If the atoms are left alone in free space in the |12; 0) state then the state will decay down 
to, 

P/.ee.pace=Tr^^|ll;0)(ll;0| + -^P— 113;0)(13;0| (27) 

i 21 i- i 23 i 21 + i 23 

which is an unentangled mixed state of the system. In fact, in free space the dressed state 
picture offers no advantage since the Rabi splitting is so small that the line widths of the 
dressed states overlap, and consequently the dressed states don't decay independently of 
each other. However, when the two atoms are placed in a cavity in which the Rabi splitting 
is larger than the line widths, the dressed states do not overlap. This gives rise to the dressed 
states \xo) and that do not couple to any of the other dressed states but themselves. 
Note that \xo) and j^^) both have the same concurrence, C — 1, and that the reduced 
density matrix of both these states are the same, namely, 

TTf,,ia{\Xo){Xo\) = Trf,M:)m) = \m) - |21»((12| - (21) ^ \r){r\- (28) 

The idea of starting out in the |12;0) state is to maintain the population in |Xo)(|0o)) 
component of the density matrix, and pumping the rest of the population out of the |x±) 
and 10") states. Why do we want to pump the population out of the \x±) states when it 
has a (7 = I? The entanglement content of \x±) comes from the Bell state compo- 
nent. Unfortunately, the incoherent mixture of different Bell state components degrade the 
entanglement, and in the worst case make the concurrence zero. The extreme case is, 

Patoms = + \r){r\) = ^(ii2)(i2| + |2i)(2i|) (29) 

which is clearly a separable state. Therefore, to preserve the entanglement in the |Xo)(|0o)) 
state we want to minimize the population in the \x±) and |(^!^) states which contain a 
component. The system we put forward does indeed pump the population out of the \x±) 
and 10" ) states and puts the population in the |33; n) state. Once the two atoms prepared in 
the 1 12; 0) state is put in the cavity the rate equations governing the population in the \x±) 
and 10") are the same as the ones we have derived in the previous section except they will 
contain a decay term that would put the population in the |33; n) state. Since the pumping 
field is not in resonance with the |3) ^ |2) transition the population in the |33; n) state will 
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not get pumped out. In the steady-state, the reduced density matrix of the atoms would 
then be, 

Patom3 = ^|33)(33| + ^|V'-)(^-|. (30) 

This is a three level state, so we cannot use Wootters' concurrence to calculate the 
entanglement content, but a simple argument can convince us that this state is entangled 
with a Bell state content of one half. Let's consider a very crude distillation protocol in 
which we project Patoms onto the |33) state. Clearly, half the time we will detect the state to 
be in |33), but the other half of the time we will get a null result in our detector suggesting 
that the atoms are in state. This means that we can obtain a Bell state half of 
the time, and therefore the state contains at least a Bell state content of one half. In this 
case, we can go a step further and say because the mixture has half its population in the 
Bell state, one cannot distill more than a half a Bell state for this density matrix. Therefore, 
our distillation protocol is indeed optimum, and there is a Bell state content of one half in 
this state. Another point to note here is that the final state is independent of the pumping, 
leakage, and spontaneous decay parameters in the rate equations. These parameters only 
determine the time scale in which the atoms reach the steady-state, but not the final state. 

Now we can ask ourselves what would happen if we started out in the |12;0) state for 
the two closed two-level atoms we have considered in the previous section. The steady-state 
reduced density matrix of the two closed two- level atoms in the cavity would be, 

Patoms = ^ IV'") + ^ {PgPg + PslPsl + Ps2ps2 + Po' ,2po' ,2) ■ (31) 

The concurrence of this density matrix is given by, 

Catoms = max{^ - -^{Psi + Ps2) - \^ {Ps2 + 2P„,,2) (2Po',2 + Ps2 + 4P, + 2P,i), O) . (32) 

This suggests that the concurrence of the closed two- level atoms can indeed be non-zero, 
and the concurrence is maximum (C = |) when all the population that is not in the |(^") or 
\Xo) is in the |11; 0)(11; 0| state. Therefore, the concurrence in this case is maximum when 
the cavity is not pumped at all. As we have noted before, the Bell state component that 
is generated through the pumping process is the l'^"'") state, and the incoherent mixture 
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of 1'^''") and l'^") degrades the entanglement between the atoms. Therefore, by pumping 
into the state the concurrence is reduced. This is consistent with the above expression 
for the concurrence between the two atoms. It follows that the entanglement generated 
between the atoms in both the open and closed two-level system is due to the asymmetric 
initial state and the presence of the cavity. The pumping itself does not directly contribute 
to the entanglement. What the pumping does in the open two-level system is that the 
pumping field puts half the population in the |33;n) state of the system to maximize the 
entanglement between the atoms. The advantage of the open two-level atoms is that if one 
were to perform operations on the two-atoms in the (|1), |2)) subspace, then the |3) state 
would not interfere with the operation while in the case of the close two-lcvcl atoms one could 
not easily determine if the operation was performed on the entangled or the unentangled 
component of the density matrix. 



IV. CONCLUSION 



We have shown that in our model for two atoms inside a cavity a pumping field cannot en- 
tangle two closed two-level atoms in the steady-state, but it is possible to have entanglement 
in the presence a pumping field for two open two-level atoms with a Bell state content of 
one half through optically pumping atomic population out of the decoherence-free subspace 
of the system. We assume a situation in which the vacuum Rabi splitting is larger than any 
of the decay parameters in the model, and we do not consider dipole-diploe interactions or 
motional effects in our calculation. The large Rabi splitting suggests the atoms and field can 
build up coherence before any decay process takes over. Our assumption allows us to treat 
the incoherent process of population transfer between the dressed state of the system with 
rate equations and account for the coherence between the atoms through the dressed states. 
This approach explicitly points out the states that contribute to the entanglement of the 
atoms (the n — 1 dressed states and the anti-symmetric dressed states). Qualitatively, one 
can see that if there are m atoms in the cavity then the dressed states associated with n < m 
excitations are the states that contain some entanglement between the atoms. Because of 
the various types of non-equivalent entanglement for m > 2 partition systems, there are 
no standard measures for entanglement and the analysis is beyond the scope of this paper. 
However, if one is only looking for a particular type of entanglement present in the system. 
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then we believe that the dressed state picture we present offers a more direct way to see 
where the target states he in the dressed state ladder. 
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